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. Abstract 

We study the representations of tensor random fields on the sphere basing on the theory 
of representations of the rotation group. Introducing specific components of a tensor field 
and imposing the conditions of weak isotropy and mean square continuity, we derive their 
spectral decompositions in terms of generalized spherical functions. The properties of 
random coefficients of the decompositions are characterized, including such an important 
question as conditions of Gaussianity. 
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^ ! 1 Introduction 

This paper was inspired and motivated by recent advances in the study of cosmic microwave 
background (CMB) radiation, which is currently a crucial topic of investigation in cosmology. 
Amazingly, the CMB anisotropies contain a wealth of important cosmological information and 



Q\ [ provide the main testing ground for the theories of early universe, with the possibility of precise 
determination of fundamental cosmological parameters. Note that this area has raised a number 
of issues, theoretical and practical, both in physics and mathematics. Nowadays, with the flood 
^ ■ of high quality data available from several satellite missions and expected from the future ones, 
the focus of CMB research is on developing theoretical tools for data interpretation and analysis. 

The CMB radiation is characterized by four Stokes parameters [9]: the intensity / (or, 
equivalently, the temperature T) and the other three parameters Q, U and V, which define the 
polarization state, namely, linear and circular polarization (however, the fourth parameter V 
describing the circular polarization is not necessary in the standard cosmological models). 

Much more attention has been paid to the study of temperature fluctuations in cosmic 
microwave background, great deal of experimental activity has been accompanied by thorough 
numerical and analytic work. In particular, rigorous mathematical basis has been provided 
through the excellent series of papers pQj, [2], [15] -[2T]. 
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From the mathematical point of view, the temperature fluctuations are considered as (a 
realization of) a random field on the unit sphere. One natural tool for statistical analysis of 
random fields is a spectral decomposition, which for scalar random fields on the sphere is a 
decomposition in the series of spherical harmonics (see, e.g., [13], [32]). These decompositions 
have been known and used in different applied areas for quite a long time, however recent 
contributions in p], [2], [16]-[2T] shed more light on many important practical issues such as 
characterization of spherical harmonics coefficients and their asymptotic behavior, character- 
ization of spectrum and bispectrum, testing for Gaussianity on spherical random fields and 
some other. Higher order spectra for spherical random fields have been analyzed in [20] . 

It has been long recognized that polarization in CMB provides a source of information on 
cosmological parameters complementary to temperature fluctuations. This significant addi- 
tional information is necessary to create a complete picture of anisotropies used to constrain 
cosmological models, to improve the accuracy in determining cosmological parameters, and 
also to study relic gravitational waves, which is one of the biggest challenges today (see, e.g., 
[34J, [3j). However, polarization anisotropies in CMB are detected at much smaller scales 
than temperature fluctuations and this limitation has been overcome only recently with im- 
proved sensitivity and resolution of CMB experiments. First polarization measurements were 
obtained by WMAP experiment (see, e.g., [25], http://map.gsfc.nasa.gov/) and further data 
are expected from Plank surveyor satellite (see http://www.sciops.esa.int/) prompting many 
groups of scientists to focus on polarization research. Last few years of intensive studies in this 
rapidly expanding area resulted in developing physical and mathematical formalism to describe, 
characterize and analyze the polarization counterpart in CMB. 

From the mathematical point of view, considering polarization we deal again with random 
fields on the sphere, however, polarization Stokes parameters Q and U are not scalars but rather 
the components of a symmetric trace-free rank 2 tensor, or, equivalently, the combinations 
Q ± ill are spin ±2 quantities. Therefore, for their mathematical description one needs to 
consider tensor-valued random fields (or random fields with spin-weighted quantities) on the 
sphere and their spectral decomposition and statistical analysis. 

In order to study polarization anisotropies, analytic methods have been formulated in phys- 
ical literature based on expansions of tensor-valued functions on the sphere in the series of 
appropriate tensor spherical harmonics, or expansions in spin-weighted spherical harmonics 
(this was done initially by two group of scientists in [12] and [33] respectively, and numerous 
further developments have been elaborated by now). The main reference for the mathematical 
formalism used in these studies comes from non-probabilistic literature, in particular, the tools 
for the mentioned above expansions are provided by the group representation theory, harmonic 
analysis on the group of rotations and the Peter- Weyl theorem. However, more caution and 
accuracy are needed when considering expansions for random functions, in particular, one needs 
to precise in what sense and under what conditions the corresponding series with random coef- 
ficients converge. These important issues are missing in physical literature, or at least they are 
not rigorously and clearly pronounced. It is also important to reveal in detail the nature and 
features of the random coefficients of the series expansions. Note that rigorous probabilistic 
framework is necessary to proceed with statistical inference methods. 

Only recently the rigorous mathematical and probabilistic basis to underpin used in physical 
literature techniques appeared in several papers: [7], [8], [15]. In these papers the theory for 
random sections of fiber bundles over the sphere is developed with cosmological applications in 
mind. 

In this paper we present an approach to derive spectral decompositions of tensor-valued 
random fields on the sphere from more group-theoretical point of view. We introduce specific 
components of a tensor field and under the conditions of weak isotropy and mean square 
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continuity define the representation of the rotation group by rotational shift operators acting 
on the Hilbert spaces generated by these components. Then spectral decompositions of the 
components are obtained as the sum of projections on the irreducible representation spaces, 
which expressed in terms of generalized spherical functions. We characterize the properties 
of random coefficients of the expansions obtained, including such an important question as 
conditions of Gaussianity. 

Note that problems involving spectral analysis for random fields on the sphere arise in 
various applied areas such as geodesy, geophysics, planetary sciences, astronomy, cosmology, 
medical imaging etc. 

The organization of the paper is as follows. We start in Section 2 with the brief overview 
of the representation theory of the rotation group, to that extend which will be needed in the 
further exposition. In Section 3 we review different approaches used to derive the spectral 
representation for scalar random fields on the sphere, with the particular attention paid to the 
approach based on representation of the rotation group with representation space, which is the 
Hilbert space generated by the random field. This last approach will be applied in Section 4 to 
derive spectral representations of vector and tensor fields. 

2 Elements of the theory of the rotation group 50(3) and 
its representations 

We summarize briefly some important facts on the group SO (3), which will be used in the 
following sections. In our exposition we prefer to follow to the classical books [6], [31], we will 
also indicate the connection to some other definitions and notions used in the literature, in 
particular, we will refer to yet another classical textbook [29]. However, many other excellent 
sources on the topic are available. 

The group 5*0(3) of special orthogonal transformations in R 3 consists of all rotations in 
R 3 about a fixed origin. This group has a convenient realization as the group of all 3 x 3 real 
matrices A such that A 1 A = I and detA = 1. 

A number of parametric representations of rotations can be introduced. One important 
parametric form for the group of rotation uses the so-called Euler angles. Note that there exist 
different conventions how to define these angles. One of such conventions (e.g., [6]) is to define 
a rotation g = g((pi,9,ip 2 ) as the product of three successive rotations: a rotation g Vl by an 
angle ipi around the z axis, then a rotation go by an angle 9 around the new x axis, then a 
rotation g V2 by an angle ip 2 around the new z axis. Correspondingly, the elements of the matrix 
of the rotation can be expressed explicitly in terms of the Euler angles (<pi, 9, <p 2 ) (see, e.g., [5]). 

Let V be a Hilbert space. A representation of a group G with a representation space V is a 
homomorphism T : g —>■ T(g) of G into the space of bounded linear operators on V. Thus, the 
mapping T satisfies the conditions: 

T( gi )T(g 2 )=T( gi g 2 ), T(e) = E, 

where e is an identity element of G, E is an identity operator. Note, that with each operator 
T(g) one can associate its matrix {tij(g)} (defined with respect to an orthonormal basis in V). 

A representation T is said to be irreducible if there does not exist a proper subspace W of 
V, which is invariant under T. Otherwise, T is reducible. The representation T is said to be 
unitary if the operators T(g), g G G, are unitary with respect to the scalar product defined on 
V, i.e., for all x, y G V, g G G we have 

(T(g)x,T(g)y) = (x,y). 
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We formulate now the basic results on the resolution of a representation into irreducible 
representations (see, [6]). 



Proposition 1 Let T : g — > T(g) be an unitary representation of the group of rotations in a 
(separable) Hilbert space V . Then there exist mutually orthogonal finite-dimensional subspaces 
Vi, V2, . . ., invariant with respect to T, in each of which the representation T is irreducible, and 
the space V is the orthogonal sum of these subspaces Vi. This means that every x G V can be 
expresses as a convergent series x = Y2 x i> x * ^ Vi> where convergence is meant with respect to 
the norm generated by the scalar product defined on V. 

Proposition 2 Each irreducible representation of the group of rotations is defined by a number 
I, called the weight of representation, the corresponding invariant representation space Di has 
dimension 21 + 1. The matrix of the representation (of weight I) corresponding to an arbitrary 
rotation g = g(tpi, 6, ip 2 ) has, in the canonical basis, the form 



T(g) = {T l mn (^0^2)}m,r. 



where the elements are given by the following functions 

T l mn {^e^ 2 ) = e-^P l mn { C o*6)e- m ^. (2.1) 

Note that the matrix representation (12.11) is referred to the fixed canonical basis in the 
space Di, which is formed by 21 + 1 orthogonal vectors (of dimension 21 + 1), each being the 
eigenvector with the eigenvalue e~ inip , n — —I, . . . ,1, for the rotation around z axis by an angle 
if. 

The matrix elements T^ n (<^i, 9, <^ 2 ) are called the generalized spherical functions of the order 

KM)- 

The functions P^ni?) appearing in the formula (12.11) can be represented in different ways. 
We give here their differential representation, which is also called the Rodrigues formuls0: 



P l (z 

mn\ 



(-ly-^m-n r (l + n )\ ] 1/2 

(2.2J 



2 l l(l-m)\(l + m)\(l-n)\ 
x (1 - z)- {n - m )/2 (l + z )-^ +m )l 2 ^- [(1 - z ) l ' m (l - z) l+m ] . 

Expressions for P^ n (z) can be given in terms of hypergeometric functions, trigonometric 
functions, there exists also their integral representation. All this can be found, for example, 
in [31], [29], where the connection of P l mn {z) with classical orthogonal polynomials and some 
other their properties are also reported. 

In many cases one considers representations of a group G by shift operators in the linear 
space of functions defined on some space X, in particular, X can be equal to G itself. These 
shift operators acting on functions f{g), g G G, can be defined as 

T{g )f{g) = figa'g), 

which corresponds to the case when G is considered as the group of left shifts, or one can define 
the action of the operator T(g ) as 

T (9o)f(g) = f(99o), 



1 The formula f|2 . 2|) is actually called the Rodrigues formula in [31] when P^„(z) is defined via (|2.2|) with 
^m-n dropped^ instead % m ~ n is included in (|2.ip . 
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in the case when G is considered as the group of right shifts. 

Let L 2 (G) be the Hilbert space of all functions on the compact group G, which are square 
integrable with respect to the Haar measure dg on G (this measure is invariant for compact 
groups). Considering the case G = SO (3), we have the space of functions 

f(9) = f(<Pi,0,<P2), 

for which the following integral exists 

Jjf(g)\ 2 dg^-^J J J\f( Vl ,9, V2 )\ 2 sm9d Vl d9d^ 2 <oo; (2.3) 
the scalar product of fi(g) and /2G7) is defied as 

(/i,/ 2 ) = J h(g)Mg)dg. (2.4) 

The transformation 

T(g )f(g) = f(ggo) 

forms a unitary representation in the space L 2 (G) (called the (right) regular representation of 
the rotation group [6]). The irreducible representations into which it can be resolved are the 
representations in the subspaces of generalized spherical functions T] nn ((p 1 , 9, cp 2 ) for a fixed I 
and m. Thus, the resolution of this representation into irreducible representations means that 
every function / e L 2 (SO(3)) can be expanded as a series in the functions T4 n (<^i, 9, ip 2 ). 

Proposition 3 The set of the generalized spherical functions T^ n (y?i, 9, (p 2 ) (I being an integer) 
forms a complete orthogonal system in the space of functions L 2 (SO(3)) . 

Note that the above statement is a particular case of the Peter- Weyl theorem, which is one 
of the most important results of the harmonic analysis on compact groups (see, e.g., [31]). We 
will return to this theorem later in this section and consider its stochastic version. 

The multiplicative law of the group representation implies the following rule according to 
which the generalized spherical functions are added. 

Proposition 4 (Addition formula for the generalized spherical functions.) 

1 

T l mn {gm) = E T L{gi)T l s M- (2.5) 

s=-l 

Using the unitarity of the representation matrices, the above formula can also be written in 
the following form 

TUg^ 1 ) = E TUgJUM- (2.6) 

s=-l 

Let us return to the definition of the Euler angles (ip 1 ,9,ip 2 ). Another way to introduce 
them, most commonly used in physical literature, e.g., in quantum mechanics, is based on 
so-called zyz-convention about rotations (the introduced above is zxz-convention), when a 
rotation is defined as to be produced in the following three steps: a rotation g^ by an angle 
(fx around the z axis, then a rotation gg by an angle 9 around the new y axis, then a rotation 
g V2 by an angle ip 2 around the new z axis. 
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Correspondingly, within this approach, the matrix elements of the irreducible unitary repre- 
sentation of the weight I of the rotation group 50(3) are represented by the Wigner D-functions 
D l mn ((pi,9,(p 2 ) (see, [30], [I], |29j). Proposition 3 holds for the functions D l mn , that is, every 
function from L 2 (SO(3)) can be expanded in a series in terms of the Wigner D-functions. 

Note that in quantum mechanics these functions appear to be the elements of the matrix 
representation of the rotation operator in the basis formed by the angular momentum eigenvec- 
tors. Wigner D-functions play an important role in various fields of modern physics including 
nuclear and molecular physics. 

Both conventions about the definition of the Euler angles are closely related and the same 
rotation R can be achieved with a simple adjustment of angles: R zyz (a, /3,j) = R zxz {a + 
7T/2,/?, 7 -7r/2) (see, |22|). 

Correspondingly, there exists a simple relation between the functions T l mn and D l mn (see, 
[29]): 

D l mn {^9^ 2 ) = {-i) n - m T l mn {^0^ 2 ). (2.7) 

However, an advantage of the Wigner D-functions is that their factorized form is given in 
terms of the real Wigner <i-functions d l mn (8) and the complex exponentials: 

D l mn {^9^ 2 ) = e-^cC(cosfl)e-^, (2.8) 

where d l mn are equal to P l mn with the factor (—i) n ~ m dropped (cf. (12. ip - (12.21) ). 

As we have seen above the group representation theory (outlined here for a particular case 
of the group G = SO (3)) provides a tool - the Peter- Weyl theorem - for decomposition of 
L 2 (G) space (endowed with the Haar measure) into an orthogonal sum of finite dimensional 
spaces and, correspondingly, the decomposition of functions from L 2 (G) into the sum of their 
projections on these spaces which can be represented in terms of appropriate basis functions. 

Our main interest is in extension of this result which allows for decomposition of random 
functions. One can observe immediately that such a decomposition can be derived for random 
functions X(g), g G G, whose trajectories g — > X(g) are P-a.s. square integrable w.r.t. the 
Haar measure (that is, belong to L 2 (G)), corresponding decomposition must be understood in 
the L2(G)-sense, P-a.s. However, more profound results can be obtained with the assumption 
of isotropy, that is, invariance in law of a random function X(y), y &Y, under the action of a 
group G (with G being a topological compact group acting on 7, F can be equal to G itself). 
In such a case one can formulate a stochastic version of the Peter- Weyl theorem (see, [26]) and 
the decomposition of L 2 (G) translates easily into construction of spectral representations for 
isotropic random functions. To be more precise, for the case of the rotation group SO(3) the 
following result can be stated as a consequence of the stochastic Peter- Weyl theorem (see, [20J , 

Proposition 5 Let X(g), g G G = SO(3), be a square integrable (L 2 (P(cLo)) , strictly isotropic 
random field, that is, finite- dimensional distributions of 

{X(g 1 ),...,X(g k }} and {X{gg 1 ), . . . , X(gg k )} 

are the same \/k, Wg, gx, . . . , g^ G G. Then 



X(g) = X{<p,0,4,) = J2J2 a ^y^T Dl mn(^0^), (2.9) 

l m,n 



6 



both in L 2 (P(du) x dg) and pointwise in L 2 (P(du)) , where dg is the Haar measure on SO (3), 
and 



aimn = J G X(g)^I^^D l mn (g)dg. (2.10) 



Remark 1 The functions D l mn (g) (as well as T l mn (g)) are orthogonal but not orthonormal. 
The orthogonality relation for D l mn (g) reads 



f &71 2 

J G D l m n(g)Di, n ,(g)dg = — i B <L.'* 



This explains the appearance of the factor y^jr in (I2.9l) - fl2.10l) . 



Remark 2 In what follows we will refer to the above result, however we will also consider the 
representation of 5*0(3) via (rotational) shift operators defined on the Hilbert space generated 
by a random field, which gives rise to alternative derivation of the spectral representation (12.91) 
and allows to treat the coefficients from the different point of view. 



3 Spectral representation of scalar fields 

In this section we review different approaches used to derive the spectral representation for 
scalar random fields on the sphere. 

In what follows we will suppose that we are given the probability space (Q, J 7 , P). 

Let L 2 (P(duj)) be the Hilbert space of random variables Y such that i?|F| 2 < oo, with 
inner product (Yi,Y 2 ) = E{Y{Y 2 )] a point in the unit sphere S 2 be denoted t = (l,9,(p); 
X(t) = X(t,u) be mean square continuous zero mean random field on 5*2, L 2 x (P(duj)) be the 
Hilbert space generated by X(t), G be the group of rotations, G = SO (3). 

1. First approach (which we present here following to [11] and [21]) is based on the theory 
of representations of the rotation group. The representation of the rotation group is defined by 
(rotational) shift transformations on the Hilbert space generated by the random field X(t). 

X(t), t e S 2 , can be regarded as a random field on G: 

7T 

X(t) = X(g t e ), t = g t e , g t E G, g t = g(<p + -, 9, ip 2 ), (3.1) 

where eo denotes a unit vector along the polar axis. (We suppose we have fixed the Cartesian 
coordinate system with the unit vectors (e x ,e y ,e z ), with e z = e being the unit vector along 
the polar axis of the sphere S2 and the origin is at the center of the sphere.) 

The scalar field (13. If) is independent of the third Euler angle or of the rotation around t. 

The field X(t) is said to be isotropic in the wide sense, if the correlation function R(ti,t 2 ) = 
E{X{ti)X{t 2 )) is invariant under arbitrary rotation: 

R(t 1 ,t 2 )=R(gt 1 ,gt 2 ),geG. (3.2) 

This implies that R(ti,t 2 ) is a function only of cos9, where 9 is the angle between the unit 
vectors t\ and t 2 with spherical coordinates (6*i,y?i) and (9 2 ,ip 2 ), 
cos 9 = cos 6*i cos 9 2 + sin 9\ sin 9 2 cos^x — ip 2 ). 
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Define the rotational transformation of the field X(t): 

S 9 :X(t)^X(g-H), geG. (3.3) 

This transformation induces a transformation on random variables from L 2 X : 

U 9 : Y -> U 9 Y, geG, Y G L 2 X . (3.4) 

From invariance (I3.2p it follows that the scalar product in L\, that is the covariance of Y, Z G 
is invariant under f/ 9 : 



E{U 9 YU9Z) = EiYZ), (3.5) 
which implies that U 9 is a unitary operator on L\. Also, we have the group properties for U 9 : 

U 91 U 92 = U 9192 , (U 9 )- 1 = U 9 ~\ U e = 1. (3.6) 

This means that U 9 , g G G, gives the unitary representation of the rotation group G in the 
Hilbert space L x . From m.s. continuity of X(t) it follows that U 9 is continuous w.r.t. g. 
In particular, U 9 X(t) = X(g -1 t), therefore, we can write: 

X(t)=X(ge ) = U 9 - 1 X(e ). (3.7) 

Denote by H the subgroup of rotations around the polar axis e and by H t the subgroup of 
rotations around the vector t. Then for the scalar filed X(t) we have 

U h X(t) = X(t), h G H t , U h X(e ) = X(e ), heH. (3.8) 

By the representation theory of the rotation group, the representation space L\ for U 9 can 
be decomposed into the sum of the irreducible spaces. Correspondingly, the vector X(e ) from 
L\ can be decomposed into the sum of the vectors of this orthogonal irreducible spaces. Denote 
by Di(Q) an irreducible space of the weight I representation for U 9 . Denote the canonical basis 
for D t (p.) by 

Zr, m = -l,...,l; (3.9) 

the orthogonality relations are: 

E(ZJ n zf)=5 w 5 mm/ . (3.10) 

By (13. 8p . X(e ) has only 0-th canonical component for each /, that is, its irreducible decom- 
position can be written only in terms of Z®\ 

oo 

X(e ) = TTF l Z?, (3.11) 

1=0 

where the series converges in L 2 (P(duj)), the expansion coefficients can be given: 

F l = E{X{e )Zf), 1 = 0,1,2,... (3.12) 

To obtain X(t) we apply U 9 1 : 

oo 

X(t) = ^F l U 9 ~ 1 Zl (3.13) 

1=0 

The final step is to represent U 9 Zf in terms of the canonical basis Z™. We know the form 
of the matrix of the representation in the canonical basis, therefore, we can write: 



u^z? = J2 TL(g- l )zr = y E Y r( e Mzr, (3.14) 

m=—l m=—l 
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substituting this into (13.131) we obtain the spectral representation of the field X(t): 



oo / 



4) = EE F^iQMK (3-15) 



=0 m=-l 



in L 2 {P{duj)) 1 here we have denoted Fi = Fi^J 
For the covariance function we obtain: 

oo I 

r{6) = E{x{t x )x(^)) = Y,\ p i\ 2 E yr(Oi^i)Yr(02,^) 

1=0 m=-l 

oo 

= ^E( 2/ + 1 )i F 'i 2p '( cos ^' ( 3 - 16 ) 

1=0 

applying the addition formula for spherical harmonics Y l m (9, tp) 

i 

PI{C0S9) = 27TT ^ Yr(^Vi)Yr(0 2 , V2 ), (3.17) 

m=— I 

where Pi is the Legendre polynomial. 

Remark 3 It is interesting to note that within the described above constructive approach for 
derivation of the spectral representation we obtain the spherical harmonics coefficients in the 
factorized form F\Z™ (see, ( I3.15P ) with random variables Z z m being uncorrelated (orthonormal) 
in a way ( 13.101) and the factor Fi is common for all Z z m , m = -I, Note that \Fi\ 2 is called 
the angular power spectrum and FiZ™ is called the random spectrum in [21] (or multipole 
coefficients, in physical literature). 

2. The spectral representation of a random field on a sphere can be obtained as a particular 
case of the stochastic Peter- Weyl theorem (see, [2D] , [26] and the end of Section 2). Indeed, 
Proposition [5] gives the spectral decomposition for a square integrable isotropic random field 
X(g) on the group 5*0(3). If we consider the restriction of X(g) on the quotient space S 2 = 
SO (3) /SO (2), then in the representation (12.91) the inner double sum will reduce to the single 

sum over m = —I, /, the functions D l mn (Lp, 9, ip) will reduce to D l m0 (Lp, 9, ip) = yj ■^-^Y( n {9, if), 

and the coefficients ai mn given by the formula (I2.10p will simplify to 



a i<. 

with ai m given by 



0, n = 



's 2 

that is, the resulting representation appears in the form 



a lm = I X{9,(f)Y l m {9,(f) sin 9 d<f d9, (3.18) 



oo £ 



x(9,^) = j2J2 ^xrid^), (3.19) 



«=0 m=-l 

where the random coefficients are given by (13.181) . 
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Remark 4 Comparing (I3.19P and (13. 15)) we see ai m = FiZ™, that is we can treat the random 
coefficients in the spectral decomposition from different points of view. Note also that E\ai m \ 2 = 
Ff is angular power spectrum, which is more commonly denoted as C\. 

3. We would like to recall that the very elegant derivation of the spectral representation 
of a random field on the sphere can be given via application of classical results from analysis 
and probability theory, just in three steps: (i) from the Funk-Hecke theorem one can enjoy 
a complete set of eigenvalues and orthonormal eigenfunctions for the covariance function of a 
mean square continuous homogeneous isotropic random field on the sphere, hence (ii) Mercer 
theorem allows to write down the decomposition of the covariance function and finally (iii) 
Karhunen theorem implies the spectral decomposition of the field itself (for more detail see, 



e.g., m, [32]). 

4 Spectral representations of vector and tensor fields 



Let now X(t), t G S 2 , be a vector random field, that is, given a fixed coordinate system, the 
filed X(t) at each point t G S 2 can be represented by its coordinates (Xi(t) , X 2 (t) , X 3 (t)) , which 
will change according to the certain rules, when the coordinate system is changed. We will also 
suppose that the field has zero mean. 

We want to consider the field X(t), which is a vector function of two variables (9, (p) G S 2 , 
as being defined on the group of rotations G, moreover we want to represent the field by 
a system comprising three functions of three Euler angles (ipi, 9, (p 2 ), functions of rotations, 
which transform into themselves under rotations. Thus, rotational shift transformations of 
each of such functions into itself will give rise to a representation of the group SO (3) and 
analogously to Section 3 will lead to corresponding spectral decompositions. To obtain such 
functions we will use ideas from [6]. 

Fix the Cartesian coordinate system with the unit vectors (e x ,e y ,e z ), with e z = e$ being 
the unit vector along the polar axis of the sphere S 2 and the origin at the center of the sphere. 
Let (e^e^e^) be the basis unit vectors for the polar (spherical) coordinate system, this triad 
is associated with the point (l,9,<p) of the unit sphere and constitutes the so-called 'moving 
frame'. With each rotation we can associate, however, the triad related to an element g of 
the rotation group. We associate with the point Pq, the 'North Pole' of the sphere, the triad 
(e x , e y , e z ) and this will correspond to the identity rotation e. Each rotation g transforms the 
triad at P into a triad (ei,e 2 ,e 3 ) and places it at some point P on the sphere, that is, the 
element g of the group 50(3) determines the triad and the point which it is associated with. 
The triad corresponding to the rotation g = g(<fi, 9, 0) is (— e v , eg, e r ) and it is placed at point 
P((p,9), where <p = | — ipi, and to the rotation g = g(ip 1 ,9,ip 2 ) there correspons the triad 
(ei,e 2 ,e 3 ) given by 



Now considering the vector X(P) at point P we resolve it with respect to (ei, e 2 , e 3 ) to obtain 
the components (Xi(g), X 2 (g), X 3 (g)). The component 



is normal to the surface of the sphere at point P(9, ip) and does not depend on the third Euler 
angle (p 2 , from the other two components (lying at the tangent plane at point P) we form two 



e 3 = e r . 



ei = -e^cosy^ - e d siinp 2 , 
e 2 = — e<p sin cp 2 + eg cos ip 2 , 



(4.1) 
(4.2) 
(4.3) 



X 3 (g) = X r (cpi, 9, ip 2 ) = X r (9, (p), (p = (f 2 - 7r/2, 



(4.4) 
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complex components 



X+(g) = X + (<p li e i <p 2 )=X 1 (g)+iX 2 (g) i 
X_(g) = X_(<p 1 ,e, ( p 2 )=X 1 {g)+iX 2 (g). 

These complex components can be represented as 

x+{g) = e^[-x^e^) + ix e {e^)i (4.5) 

X4g) = e-^{-X v (6^)-iX e (6,v)}, (4.6) 

where X (/3 (6 I , if), Xg(8, ip) are the components of the vector field in the polar coordinate system, 
(p=*-<pi. 

With the above three functions we have defined the vector field on the rotation group: 

X(<p 1 ,0,<p 2 ) = X(g) = (X + {g),X4g),X r (g)). 

With each rotation go the normal component X r (g) = X r (P) at point P((p,6) transforms into 
the normal component at the point g^ g = g^ P and does not depend on the third Euler 
angle, that is, on the rotation in the tangent plane. 

We can consider the normal component ( 14.4D and the component of the field in the tangent 
plane, which is represented by two functions (I4.5p - fl4.6l) . separately. 

Remark 5 Note that in applications vector fields on the sphere are usually considered as 
vectors in the tangent plane at each point of the sphere. Typical example of vector fields on 
the sphere that arise in practice are electromagnetic fields or wind velocity. 

The notion of isotropy (in weak or second-order sense) we will define componentwise, as 
isotropy property for the components fl4.4p - fl4.6l) . that is, the covariance functions 

R+(gi,g 2 ) = EX + ( gi ) X + (g 2 i (4.7) 
R-( gi ,g2) = EX4 9l )X_(g 2 ), (4.8) 
R r ( gi ,g 2 ) = EX r ( 9l )X r (g 2 ) (4.9) 

are invariant with respect to rotations g G 50(3): 

Ri(gi,g 2 ) = Ri(ggi,gg 2 ), Vg, i = ±,r, (4.10) 

which implies R i {g 1 ,g 2 ) = R i (g^ 1 g 2 ), i = ±,r. 

Remark 6 We believe that defining the notion of isotropy in such a specific way, that is, 
isotropy for the components taken in the moving basis (14. ip ~ 04.31) is quite reasonable for vector 
fields on the sphere. This will be quite analogous to the notion of isotropic vector fields 
introduced in early works on turbulence (see, e.g., [27], [22] and references therein). Recall that 
in those works isotropy for a vector field in strict or wide sense is defined as a property which 
prescribes that probabilistic characteristics such as probability distributions for the values of 
a field in a system of points (strict sense) or just correlations (wide sense) are invariant under 
translations (that is, motions, which include rotations, shifts and reflections) of the the system 
of points performed simultaneously with the same movement of the coordinate system. In 
particular, for the correlation tensor of the vector field this definition entails the following 
property: 

B(t) = G*B(Gt)G, 
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where G is a matrix of the corresponding transformation and G* is its transpose. (In some 
studies a field with the above property of its correlation tensor is also called isotropic in a 
vector sense.) Considering random fields on the sphere we are restricted to the rotations 
only and define the isotropy as invariance of probabilistic characteristics under rotations. In 
our approach we consider the vector field on the sphere defined at each point of the sphere 
via its components in the moving basis, which is associated with a rotation g. Therefore, 
these movements (rotations) of the coordinate system are already incorporated into such a 
representation of a vector field. And each component of the field transforms into itself under 
rotation. That is why it looks quite natural to define isotropy via the relations (14.101) . Moreover, 
this will give us the possibility to define the scalar product, which will be invariant under 
rotations, in the corresponding Hilbert spaces of random variables, and, therefore, we can 
parallel all the constructions and reasonings of the previous section, where scalar fields have 
been considered. 

First, as we noticed above, the normal component of the field X r (9, ip) does not depend on 
the third Euler angle and can be treated completely in the same manner as the scalar field, 
and the spectral representation of the form (13.151) can be written, that is, the representation 
via the spherical harmonics YJ m (#, and the corresponding representation of the correlation 
function in the form (13. . 

Consider now the tangent component of the filed given by two functions (I4.5l) - (l4.6p . To 
obtain the spectral representation we will follow the lines of Section 3. 

Let L x+ be the Hilbert space generated by X + (g). 

Define the rotational (or shift) transformation: 

S 9 : X + (g ) - X+GrV)), g,g e G. (4.11) 

This transformation induces a transformation on random variables from L x : 

U 9 : Y -> U 9 Y, g G G, Y G L 2 X+ . (4.12) 

From invariance (j4.7p it follows that the scalar product in L x+ , that is the covariance of 
Y, Z G L x+ is invariant under U 9 : 

E(U 9 YWZ) = E(YZ), (4.13) 

which implies that U 9 is a unitary operator on L x . Also, we have the group properties for 
jjg. u 9l U 92 = U 9192 , (U 9 )- 1 = U 9 ~\ U e = 1. This means that U 9 , g G G, gives the unitary 
representation of the rotation group G in the space L x . Assuming mean square continuity of 
X+(t), we obtain that U 9 is continuous w.r.t. g. 

In particular, U 9 X + (g Q ) = X + (<7 _1 <7o), therefore, we can write: 

X + (g) = X + (ge) = U^X^e), (4.14) 

where e is identity element of G, and 

X+(e) = X + (e ) = X x (e ) + iX y (e ) (4.15) 

that is, it reduces to the component of the field at the point eo = Po (the pole of our sphere). 

Denote by H eo (<p) the subgroup of rotations around the polar axis eo and by H g (<p) the 
subgroup of rotations around the vector geo- Then for the filed X + (g) we have 

U h X + (g) = e*<>X + (g), h G H g (<p), U h X + (e ) = e^X + (e ), h G # e ». (4.16) 
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The representation space L\ for U 9 can be decomposed into the sum of the irreducible 
spaces. Correspondingly, the vector X + (e) from L 2 X can be decomposed into the sum of the 
vectors of this orthogonal irreducible spaces. Denote by Di(Q) an irreducible space of the weight 
I, let 

Z+, m = -/,...,/, (4.17) 

be the canonical basis for Di(Q), that is the basis composed from the eigenvectors of U h , 
h G H eo ({p) (rotation about e ) corresponding to the eigenvalues e zipm , for these vectors the 
orthogonality relations hold: 

E{Z+$J) = 8 lv 5 mm , (4.18) 

As can be seen from (14.161) . X{cq) has only 'l'-st component non-zero in the basis (I4.17p . 
for each I, therefore, its irreducible decomposition can be written only in terms of Z^: 

oo 

X,(e )=$>+Z+ (4.19) 

1=0 

where the series converges in L 2 (P(du)), the expansion coefficients can be given: 

F+ = E(X + (e )zJ), I = 0, 1, 2, . . . (4.20) 
To obtain X(g) we apply U g : 

oo 

X + {g) = Y, F ? U ^ Z t ( 4 - 21 ) 

1=0 

Now it is left to notice that 

U 9 "z+= i2 T ^(9- 1 )ZL (4-22) 

m=— I 

substituting this into (14.211) we obtain the spectral representation of the field X+(g): 

oo I 

x +(g) = EE F i' T L{g~ 1 )zL (4-23) 

1=0 m=-l 

in L 2 {P{duj)). 

For the covariance function we obtain: 



R+(gi,92) = E(x + ( gi )x + ( 92 )) = J2\Fi + \ 2 E TUg^njg; 1 ) 

1=0 m=-l 

oo 



1=0 



applying the addition formula ( 12.61) and the orthogonality relations (14.181) . 
Analogously we obtain the spectral representation of the field X_(<?): 



oo I 



X -(9) = EE FfT-imig-^rm (4-25) 



1=0 m=-l 



in L (P(du)), where r.v. Z lm satisfy the orthogonality relations E(Z lm Z Vml ) = 5[i>5 mn 
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For the covariance function we obtain: 



R4 gi ,g 2 ) = E{X.{ gi )X^g 2 )) = Y,\K\ 2 E ^-i.mter 1 )^^ 1 ) 

1=0 m=-l 

oo 

= E 1^1^-1,-1(^2) (4.26) 



z=o 



Let us summarize the above reasonings in the following theorem. 

Theorem 1 Let X(t), t = (8,cp) G S 2 , be a zero-mean vector random field on the unit sphere 
and (X^O, <p),Xg(9, ip),X r (9, (p)) be its components in spherical coordinate system. Introducing 
a specific local coordinate system h4.1\ )- [flfy , which at every point of the sphere is related to an 
element of the group of rotations, consider the representation of the field X by the components 
defined on the group of rotation: two complex components X + (g) andX^(g) given by (|7?5p-([Jl|) 
and the normal X r (g) = X r (t). Suppose that these functions X + (g), X^(g) andX r {g) are weakly 
(second-order) isotropic and mean square continuous. Then the functions X + (g) and X_(g) 
can be expanded in the series of generalized spherical harmonics l[4-23\ ) an d l!l4-25\ ), respectively, 
with uncorrelated coefficients, convergence is meant pointwise in L 2 (P(du)). Corresponding 
covariance functions have series representations \4-2>4 ) and ^4-26 ). The normal component can 
be treated and expanded as a scalar random field. 



Let us turn now to the tensor fields on the unit sphere. We will consider the case of tensors 
of second rank. 

To find the expansions for tensor fields we can apply the reasonings similar to the above 
ones and represent the tensor fields by means of components (functions on the group 50(3)) 
which transform into itself under rotations and which will be multiplied by e mcp , m = 0,1,2, 
under rotations about the axis normal to the surface of the sphere. 

To obtain such components we note that components of a tensor of the second rank are 
transformed under a rotation in the same way as a product of the components of two vectors. 
Thus, for tensor random field on the sphere X(t), t G S 2 , we come to the following nine 
components - functions of rotations g = g(ifi, 9, <p 2 ): 



X rr , X w + Xqq ± i (X V — Xq^) , (4.27) 

(-X vr - iX er ) e^ 2 , (-X rv - iX re ) e lip \ (4.28) 

{-X w + iX 8r ) e- lLp \ {-X rv + iX r8 ) e~ ilfl2 } (4.29) 

X w -Xee + i (X e ^ - X^g) e 2 ^ 2 , (4.30) 

X w -Xgg-l(Xg v -X v g)e- 2i ^. (4.31) 



With cosmological applications in mind, we restrict our consideration to the tangential 
counterpart of the field, that is, we will be interested in the components (I4.30p - fl4.3ip only. 
Moreover, we will suppose that this tangential counterpart (the components of the tensor 
relative to the basis in the tangent plane) forms a symmetric trace-free tensor, that is, X w + 
Xgg = 0, X 8ip = X^g . In such a case the components fl4.30p - fl4.3ip reduce to the form 

2(X w + tX dv )e 2 ^ =A(jg), (4.32) 
2 (X w - iXg 9 ) e~ 2 ^ 2 = B(g). (4.33) 
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Now to derive the spectral representation for the random functions fl4.32p - fl4.33p we apply 
the same scheme as used above for the components (14 . 5 1) - (14. 6 p of a vector field. 

We suppose that the random functions A(g) and B(g) are weakly (or second-order) isotropic, 
that is, their covariance functions satisfy 

Ri(9i,92) = Ri(ggi,gg 2 ), Vg,gi, 92 e SO(3), i = A,B, (4.34) 

which implies R i (g 1 ,g 2 ) = Ri(g^ l g 2 ), i = A,B. 

Then we consider the representations of the rotation group by rotational transformations 
defined on the Hilbert spaces L\ and L 2 B generated by the random functions A(g) and B (g), 
which will be: (1) unitary due to the invariance of the scalar product under rotational trans- 
formations (see, (14.341) ). (2) continuous under the assumption of mean square continuity of our 
random field. 

We notice also that under rotations by an angle (p about the axis normal to the surface 
of the sphere (that is, around e r ) the functions A(g) and B(g) will multiply by e l2ip and e~ i2ip 
respectively. 

The representation spaces L\ and L 2 B can be decomposed into the orthogonal sums of the 
irreducible spaces of weight I, Df and Df , introducing the canonical bases {Zf^} and {Z^}, 
m = —/,...,/, in Df and Df correspondingly, we come to the representations for the functions 
A{g) and B(g). 

Namely, we obtain: 

oo I 

A (9) = E E MT l 2m {g- x )Zt, (4.35) 

1=0 m=-l 
oo I 

B (9) = E E ZiT-^g-^ZL (4-36) 

1=0 m=-l 

where the series converge in L 2 (P(du)), r.v.'s Zf^ and Zf m satisfy the orthogonality relations 



E(Z\ m Z\,J) = SwSram', i — A,B, (4.37) 
the coefficients Ai and Bi are given by the formulas 



Ai = E{A{e )Z£), Bt = E(B(e Q )Z%), I = 0, 1, 2, . . . (4.38) 
For the correlation functions we obtain 



Ra(9i 1 92) = E \M 2 T l 22 {9i l 92), (4.39) 



1=0 
oo 



R B {gx X g 2 ) = E \Bi\ 2 T l _ 2t _ 2 (g^g 2 ). (4.40) 

1=0 

We come to the following theorem. 

Theorem 2 Let X(t), t = (9, if) G S 2 , be a zero-mean tensor-valued random field on the 
unit sphere. Suppose that the rank of the tensor is 2, its components are given with respect to 
polar coordinate system and tangential counterpart of this tensor forms a symmetric trace-free 
tensor. Consider the representation of the tangential counterpart of the field X by two complex 
components defined on the group of rotation: A(g) and B(g) given by l[4-32\ )- [4.33\ ) and suppose 
that these functions are weakly (second-order) isotropic and mean square continuous. Then the 
functions A(g) and B(g) can be expanded in the series of generalized spherical functions ^4-35^ 



and Iji4.36\ ), respectively, with uncorrelated coefficients and convergence is meant pointwise in 
L 2 (P(duj)). Corresponding covariance functions have series representations ( 4-39^ and ^JD ). 
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Remark 7 We restrict ourselves here to consideration of representations for the random func- 
tions fl4.32l) - fl4.33p . Note that the same representation holds for more general functions (I4.30P - 
(I4.3ip . The functions (I4.27P can be expanded in terms of the usual spherical harmonics Y^,(0, <p), 
the functions (14.281) - in terms of T[ m , and the functions (14.291) - in terms of T l _ x m . Moreover, 
in the similar way the representation of tensor fields of higher ranks then 2 can be obtained 
(see also Remark [9] below) . 

Let us consider more carefully the random coefficients Z^, Z^ m and Zf^, Zf m of the expan- 
sions (I4.23I) - (I4.25I) and (I4.35j) - (I4.36[) . First immediate properties of these random variables we 
summarize in the following lemma. 

Lemma 1 Under the condition of mean square continuity and weak isotropy of the random 
function X + (g), the following statements hold: 

(1) the random variables Z^ m are uncorrelated with zero mean and unit variance: 

E ( Z lm Z l/m>) = $U'$mm>] 

(2) for all l,Z^= (-l) l - m Z+_ m , and, therefore, E(Z+ m Zj m ,) = for all (Z,m) ^ (l',m'); 

(3) for all l,geG, Z+ m I Y!s=-i TLWL m Particular, Z+ = e im *Z+ m V<p. 

The same properties hold for the random variables Z lm and Zf^, Zf m under the conditions 
of mean square continuity and isotropy of the corresponding fields. 

Proof. All three facts follow from the very definition of the random variables Z^ m1 which have 
been chosen as canonical bases in the spaces D\. For every /, {Z^, m = —I, ...,/} are orthonor- 
mal bases in orthogonal spaces Di, which means that (1) holds. Moreover, {Z^, m = —I, ...,/} 
constitute the canonical bases, for all I, therefore U h Zf m = e m9 Zf my h being the rotation 
around e$ by an angle ip. Taking into account that U 9 is an unitary operator for which the con- 

travariant basis jz^ m \ m = —I, / j coincides with j-Z^, m = —I, Z j and using the re- 
lation between the vectors of canonical basis and its cotravariant basis: Zi™' = (— l) l ~ m Z^_ m 
(see [H]), we come to the statement (2) of the lemma. Further, we know that for every I, 
T l (g) = {T^ s (p), m,n = —I, l\ are the matrices of the representation U 9 in the spaces Di, 
with respect to the canonical basis in Di, that is with respect to {Z^, m = —I, ...,/}. This 
implies the equality U 9 Z^ n = Yl l s =~i T} ns (g)Z^ m . On the other hand, isotropy of the field X+(g) 

entails the equality U 9 Z^ m = Z^ m VZ, Vm. This proves the statement (3) of the lemma. Since 
all the random variables Z^ m , Z^Zj^, Zf m have the same nature, the properties stated in the 
lemma are common for all of them. 

Remark 8 Considering 21 + /-dimensional vectors composed from the random variables Z^ m 
with fixed /, that is the vectors Z^ = {Z^, m = —I, we can reformulate the statement 

(3) of Lemma 1 in the following form: Zf = T l (g)Zf. Statement (1) of Lemma 1 can be 
reformulated in a slightly more general form as E(ZfZf) = hi+i and E(Z^Zjt) — 0, I I', 
where I21+1 is (2/ + 1) x (21 + 1) identity matrix and denotes (21 + 1) x (2V + 1) zero matrix. 

The random variables Z^ m , Z^ m and Z^, Zf m can be characterized further in the similar 
manner as this was done in [2] for the coefficients of expansions of scalar fields on the sphere. 
The properties stated in Lemma [1] is all that one needs to apply the same reasonings as in [2] . 
Although the results appear to be completely the same as for the scalar fields, we present them 
here for completeness of the exposition and because it is really remarkable to get the same 
characterization of random coefficients within a more general construction than scalar fields. 



16 



Lemma 2 Under the condition of mean square continuity and weak isotropy of the random 
function X + (g), the the random variables Z^ m have the following properties: ReZ^ m = ImZ^, 

the ratio is distributed accordingly to a Cauchy distribution; ReZ^ m and ImZf m are 

Im 

uncorrelated with variance \; the marginal distributions of ReZ^ m and ImZ^ m are symmetric: 
ReZ^ m = —ReZ^ m , ImZ^ m = —lmZ^ m . The same properties hold for the random variables Z~[ m 
and Z{^, Zf m under the conditions of mean square continuity and isotropy of the corresponding 
fields. 

Lemma 3 Let the random function X + (g) be mean square continuous and weakly isotropic. 
If we assume in addition that X + (g) is Gaussian, then Z^, I = 0,1,..., m = —I,. ..I, are 
independent Gaussian variables. The same is true for the random variables Z^ m and Z^, Zf 
under the corresponding conditions on the underlying fields. 



in 



Note that to prove the independence in the above lemma the statement (2) of Lemma [T] is 
essential. 

Finally, the following very important theorem can be stated. 

Theorem 3 Let the random function X + (g) be mean square continuous and weakly isotropic. 
Then for all I, the coefficients {Z^, m = are independent if and only if they are 

Gaussian. The same is true for the random variables Z^ m and Z^, Zf m under the corresponding 
conditions on the underlying fields. 

The proof is based the Skitovich-Darmois theorem, which gives the criterion for Gaussianity 
of a collection of independent random variables via independence of some linear statistics of 
these variables. For detail we refer to [2]. 

Let us return to the representations (I4.35P and (I4.36p . For the rotation g = g(<fi, 0, (p 2 ) we 
know that the inverse rotation g~ l is given by the Euler angles (n — tp 2 , 6, it — p\). Taking into 
account the expression for T l mn (see (12.1)) ) and canceling the factor e ±2tip2 in both sides of (jH 
and (I4.36P respectively, we come to the representations 

oo I 



2 (X w (6, <p) + iX 6v (e, <p)) = E AlT ^ M - <Px)Zh (4.41) 

1=0 m=-l 
oo I 



oo I 



2 (X vv (9, ip) - iXtopiO, ip)) = E ^ T -2, m (0, e, 7T - p l )Zf m (4.42) 

1=0 m=-l 
oo I 

= EE B l TL 2tm (Q,0,*/2 + <p)Z% n 



(we recall that tp — ir/2 — (pi). 

We can also reformulate the above results using Wigner D-f unctions. 

The spectral representation for the fields A(g) and B(g) can be obtained as special cases of 
Proposition O More precisely, due to the difference in the convention for rotations, we have to 
consider 

2(X ee + tX e ^e- 2i ^ = A\g), 
2{X ee -iX ev )e 2 ^ = B\g). 
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Due to the special form of A'(g) and B'(g), in the formula (12. 9p only the functions D l m _ 2 (g) 
and D l m2 (g) will participate respectively. To use Proposition [5] we need also to adjust the 
conditions imposed on A'(g) and B'(g). Namely, we demand the fields A\g) and B\g) to be 
square integrable and strictly isotropic (as described in Proposition [5]). Then we can write the 
decompositions 

A '(9) = EE <W ^r D k-2(9), (4-43) 
;=0 m=-l 

00 1 lo] 4- 1 

^) = EE b ^^r Dl mAdl ( 4 - 44 ) 

i=0 m=-Z 

where 

a lm = J A\g)^^^Di^r 2 (g)dg, (4.45) 

6 im = J B'(g) ^±lD^ 2 (g)dg, (4.46) 
the convergence of series is in Li2(P(du) x dg) and poinwise in L,2(P(duj)). 
We can also deduce from fl4.43l) - fl4.44p the representations 

oo I nrz 7" 

2 (Xee(B, tp) + iX 6lp (9, ^)) = E E a <™ V -^-^m,- 2 (^ ^ °)» ( 4 - 47 ) 
2(X et {6,<p)-iX e¥) (6,<p)) = EE ^J^^^, 2 (^^0). (4.48) 

i=0 m=— i 



Remark 9 We can see that the key step to derive the spectral decomposition for vector and 
tensor fields on the sphere is to consider the components which transform in a specific way under 
rotations about the axis e r . In harmonic analysis on the sphere there exists a special class of 
functions which are defined relatively their behavior under the rotations of the basis vectors in 
the plane tangent to the sphere (i.e. rotations around e r ). These are so-called spin- weighted 
functions. We recall the definition of these functions (see, [10], [23], |28j), however here we use 
different convention for rotation than the above authors to be consistent with the literature 
on CMB polarization. First of all, we note that these functions, similarly to the vector- or 
tensor-valued functions, are defined relatively to some coordinate system. Let at any point of 
the sphere we have defined three orthogonal vectors: one normal (radial) e r and two tangential, 
(ei,e2). A function f(9,(p) defined on the sphere S2 is said to have spin s if under a right- 
handed rotations of basis vectors (ei, e%) by an angle ip it transforms as f'{9, <p) = e~ ls ^ f(9, if). 
Spin functions can be equivalently defined as evaluation at ip 2 = of any function in L 2 (SO(3)) 
resulting from an expansion for fixed index n in the Wigner D-functions D l mn (ipi, 9, (p 2 )- Thus, 
the functions D mn ((p,9,0) or D 1 ^ __ n (<p, 9, 0) define an orthogonal basis for the expansions of 
spin n functions in i^S^). After normalization in L 2 (S2), these basis functions, which are 
called the spin- weighted spherical harmonics of spin n, are given in a factorized form in terms 
of the real Wigner c/-functions d l mn (9) and the complex exponentials e tmip as follows: 



= \/^<-n(cos^)e^. (4.49) 
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We refer for the rigorous mathematical theory for spin- weighted functions to [7], [8]. 
For example, X$g(d,(p) ± iX$tp(9,<p) are spin ±2 functions and the decompositions (14. 47ft - 
(I4.48P are decompositions in spin ±2 spherical harmonics. 

Example. In cosmological studies the observable linear polarization field is described in terms 
of two Stokes' parameters Q and U defined in a given direction with respect to the particular 
choice of axes in a plane perpendicular to the direction of observation. When reference frame 
is rotated around the direction of observation (direction of propagation) Q and U transform 
like the components of a 2-dimensional 2-d rank symmetric trace-free tensor. Therefore, linear 
polarization can be conveniently described as a tensor-valued field on a sphere 

If Q U \ 
2\U -Q ) 

defined with respect to the polar basis vectors (eg, e v ), or, equivalently, via combinations Q±iU. 

Using the above theory, we can expand Q(6, (p)±iU(6, <p) in the spin ±2 spherical harmonics. 
For this we must impose the assumption of strict isotropy, as in Proposition [5j then we can 
write 

oo I 

p(e, <p) = Q(e, if) + iu{e, <p) = E ^ ( 4 - 5 °) 

1=0 m=-l 
oo I 

p*(o, <p) ee Q(e, <p) - iu(9, <p) = E -m™ - 2Y ™& ( 4 - 51 ) 

1=0 m=-l 

both in L 2 (P(gLj) x dg) and pointwise in L 2 (P(dLj)). 

However, as our first approach suggests, the above decompositions are also hold for mean 
square continuous random fields under the assumption on weak (second-order) isotropy point- 
wise in L 2 (P(duj)). From the theory presented in this section we have some characteristic prop- 
erties of random coefficients in the above series (for mean square continuous weakly isotropic 
random fields). We know, in particular, that the coefficients are uncorrelated, for Gaussian 
fields the coefficients are Gaussian independent variables, and the coefficients are independent 
if and only if they are Gaussian. Thus, with the rigorous probabilistic framework provided in 
the paper, expansions of the form (I4.50p - (l4.5ip . which are widely used in the literature on CMB 
polarization, now become operational. 

References 

[1] Baldi, P., Marinucci, D., Varadarajan, V. S. On the characterization of isotropic Gaussian 
fields on homogeneous spaces of compact groups, Electron. Commun. Probab., 12, 291-302 
(2007). 

[2] Baldi, P., Marinucci, D. Some characterizations of the spherical harmonic coefficients for 
isotropic random fields, Stat. Probab. Lett., 77, No. 5, 490-496 (2007). 

[3] Baskaran D., Grishchuk L.P., Polnarev A.G. Imprints of relic gravitational waves in cosmic 
microwave background radiation, Phys. Rev. D, 74, No. 8, 083008 (2006). 

[4] Brink, D.M., Satchler, G.R., Angular momentum. Oxford Univ. Press (1975). 

[5] Cabella, P., Marinucci, D. Statistical Challenges in the Analysis of Cosmic Microwave 
Background Radiation, Annals of Applied Statistics, in press, arXiv:0807.1816 



19 



[6] Gelfand, I.M., Minlos, R.A., and Shapiro Z.Ya., Representations of Rotation and Lorentz 
Groups and Their Applications. Pergamon Press, Oxford (1963). 

[7] Geller, D., Marinucci, D., Spin Wavelets on the Sphere. larXiv:0811.2"935l 3 (2008). 

[8] Geller, D., Lan, X., Marinucci, D. Spin needlets spectral estimation. larXiv:0 907.33699vl 
(2009). 

[9] Chandrasekhar, S., Radiative Transfer. New York, Dover (1960). 

[10] Goldberg, J.N., Macfarlane, A.J., Newman, E.T., Rohrlich, F., Sudarshan, E.C.G., Spin-s 
spherical harmonics and d, J. Math. Phys., 8, 2155-2161 (1967). 

[11] Hannan, E., Multiple Time Series. Wiley, New York (1970). 

[12] Kamionkowski, M., Kosowsky, A., Stebbins, A., Statistics of cosmic microwave background 
polarization, Phys. Rev. D, 55, No. 12, 7368-7388 (1997). 

[13] Leonenko, N.N., Limit Theorems for Random Fields with Singular Spectrum. Kluwer 
Academic publishers, Dordrecht (1999). 

[14] Lyubarskii, G. Ya., The Application of Group Theory in Physics. Pergamon Press, New 
York (1960). 

[15] Malyarenko, A., Invariant random fileds in vector bundles and application to cosmology. 
ArXiv: 0907.4620vl (2009). 

[16] Marinucci, D., Regression methods for testing Gaussianity on a spherical random field, 
Random Oper. Stoch. Equ., 13, No. 4, 313-324 (2005). 

[17] Marinucci, D., A central limit theorem and higher order results for the angular bispectrum, 
Probab. Theory Relat. Fields 141, No. 3-4, 389-409 (2008). 

[18] Marinucci, D., High Resolution Asymptotics for the Angular Bispectrum of Spherical 
Random Fields, Ann. Stat., 34, No.l, 1-41 (2006). 

[19] Marinucci, D., Peccati, G., High frequency asymptotics for subordinated stationary fields 
on an abelian compact group, Stochastic Processes Appl., 118, No. 4, 585-613 (2008). 

[20] Marinucci, D., Peccati, G., Representations of 5*0(3) and angular polyspectra. Preprint. 
larXiv:0807.0687v l. 

[21] Marinucci, D; Piccioni, M., The empirical process on Gaussian spherical harmonics, Ann. 
Stat. 32, No. 3, 1261-1288 (2004) 

[22] Monin, A, Yaglom, A., Statistical Fluid mechanics: Mechanics of turbulence. MIT Press, 
Cambridge, 1971: Vol.1, 1975: Vol. 2. 

[23] Newman, E.T., Penrose, R., Note on the BondiMetznerSachs group, J. Math. Phys., 7, 
863-870 (1966). 

[24] Ogura, H. Representations of the Random Fields on a sphere. Mem. Fac. Eng., Kyoto 
Univ, 52, no 2, 81-105 (1990). 



20 



[25] Page, L. et al., Three- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Polarization Analysis, Astrophys. J. Suppl., 170, 335-376, (2007). 

[26] Peccati G., and Pycke, J.-R., Decompositions of stochastic processes based on irreducible 
group representations. Preprint. |Arxiv:math /0509569vl. 

[27] Robertson, H.P., The invariant theory of isotropic turbulence, Proc. Camb. Philos. Soc, 
36, 209-223 (1940). 

[28] Thorne, K.S., Multipole expansions of gravitational radiation, Rev. Mod. Phys., 52, No. 
2, 299-339. 

[29] Varshalovich, D.A., Moskalev, A.N. and Khersonskii, V.K., Quantum Theory of Angular 
Momentum. World Scientific Press, Singapoure (1988). 

[30] Wigner, E. P., Group Theory and its Application to the Quantum Mechanics of Atomic 
Spectra. Academic Press, New York (1959). 

[31] Vilenkin, N.Ja., Klimyk, A.U., Representation of Lie Groups and Special Functions. 
Kluwer, Dordrecht (1991). 

[32] Yadrenko, M.I., Spectral Theory of Random Fields. Optimization Software Inc., New York 
(1983). 

[33] Zaldarriaga, M., Seljak, U., All-sky analysis of polarization in the microwave background, 
Phys. Rev. D, 55, No. 4, 1830-1840 (1997). 

[34] Zhao W., Baskaran D., Grishchuk L.P., On the road to discovery of relic gravitational 
waves: The TE and BB correlations in the cosmic microwave background radiation, Phys. 
Rev. D, 79, No. 2, 023002 (2009). 



21 



